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We  describe.^  systems  of  conservation  laws  with  the  property  that  the  shock 

Xi 

and  rarefaction  curves  coincide.  -We-  givej.new  examples  of  such  systems.  These 

A 

systems  isolate  and  separate  many  nonlinear  aspects  of  shock  waves. 

The  author  believes  that  analytical  techniques  which  handle  problems  of 
uniqueness,  continuous  dependence  or  convergence  of  finite  difference  schemes  in 
these  systems,  would  isolate  components  in  a  corresponding  analysis  required  for 


general  systems  of  conservation  laws. 
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SIGNIFICANCE  AND  EXPLANATION 


The  Cauchy  problem  for  a  system  of  conservation  laws  in  one  dimension  is  the 
probl em 


+  f(u)x  -  0  , 

(1) 

u(x,0)  «  u0(x)  , 


where  u  «  (uj,  . ..,  un)  5  u(x,t)  ,  f  •  (f1#  ...»  fn)  and  *  e  R,  t  >  0.  Such 

equations  arise  in  applications  of  gas  dynamic s ,  elasticity,  oil  reservoir 

simulation  and  other  areas  of  engineering  when  diffusion  is  neglected.  In 

solutions  of  (1),  information  travels  at  speeds  given  by  the  eigenvalues  of 

3f 

the  matrix  i  ■  1,  n.  Since  this  matrix  depends  on  u,  the  speeds 

X^  depend  on  the  solution,  and  this  leads  to  the  formation  of  discontinuities 
called  shock  waves.  Shock  waves  introduce  many  nonlinear  phenomena.  For  this 
reason,  fundamental  issues  such  as  uniqueness  and  continuous  dependence  of 
solutions  on  the  initial  data  remain  unresolved.  One  severe  aspect  of  the 
nonl inearities  of  shock  waves  is  that  in  general ,  when  two  shock  waves  of  the 
same  family  collide,  waves  in  other  families  spontaneously  appear.  In  this  paper 
we  describe  systems  of  conservation  laws  which  arise  in  applications,  and  for 
which  this  aspect  of  the  nonlinearities  of  shock  waves  is  not  present.  New 
examples  of  such  systems  are  given.  These  degenerate  systems  isolate  and 
separate  nonlinear  effects,  the  author  believes  that  certain  analytical 


techniques  developed  for  these  systems  would  isolate  components  in  a 


corresponding  analysis  required  for  general  systems  of  conservation  laws. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 
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DEGENERATE  SYSTEMS  OF  CONSERVATION  LAWS 


Blake  Temple 

We  summarize  and  simplify  some  of  the  results  in  [17,18]  which  give  a  theory 
for  the  systems  of  conservation  laws  which  have  the  property  that  the  shock  and 
rarefaction  curves  coincide.  These  systems  have  milder  nonlinearities  than 
systems  of  conservation  laws  in  general,  and  represent  the  simplest  setting  in 
which  hyperbolic  singularities  can  appear.  Moreover,  examples  of  such  systems 
arise  in  many  applications  [1,4,5,6,8,12,14,16].  For  these  systems,  as  is  the 
case  for  a  scalar  conservation  law,  the  nonlinearities  affect  the  wave  speeds, 
but  no  reflected  wave  appears  after  interactions  [cf.  15].  The  wave  speeds  are, 
however,  coupled  in  these  systems,  a  phenomenon  not  present  in  scalar 
conservation  laws.  Thus  systems  with  coinciding  shock  and  rarefaction  curves  are 
a  reasonable  place  for  studying  some  of  the  problems  that  appear  difficult  in  the 
general  setting  of  conservation  laws. 

The  shock  and  rarefaction  curves  coincide  in  a  characteristic  family  if  the 
rarefaction  curves  are  straight  lines  (line  field)  or  level  curves  of  the 
eigenvalue  (contact  field).  Thus  for  2x2  systems  of  conservation  laws,  there 
are  generically  three  classes  of  systems  with  coinciding  shock  and  rarefaction 
curves:  those  with  two  contact  fields  (class  I),  those  with  one  line  and  one 
contact  field  (class  II)  and  those  with  two  line  fields  (class  III).  Moreover, 
under  suitable  assumptions,  these  systems  cannot  be  transformed  into  equations 
which  are  either  linear  or  uncoupled. 
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Below  we  give  conditions  on  the  flux  functions  which  are  necessary  and 
sufficient  for  determining  the  systems  in  the  above  classes  (we  record  here  the 
fact,  omitted  in  [17,18],  that  there  exist  nonlinear  coupled  systems  with  two 
contact  fields  which  are  in  conservative  form). 

In  [ 19]  it  is  proved  that  for  systems  lying  outside  the  class  of  systems 
with  coinciding  shock  and  rarefaction  curves,  the  solutions  are  not  L1- 
contractive  in  any  metric.  It  would  be  interesting  to  locate  the  systems  within 
the  above  classes  which  are  L^ -contractive. 

Continuous  dependence  of  solutions  on  the  initial  data  is  not  known  for 
general  systems  of  conservation  laws.  It  would  be  interesting  to  obtain 
continuous  dependence  or  uniqueness  results  for  systems  in  class  I,  II  or  III. 

In  [11,15]  it  is  shown  that  certain  finite  difference  schemes  are  total 
variation  diminishing  for  systems  in  class  III;  in  classes  I  and  II  more  serious 
numerical  diffusion  occurs,  and  this  result  is  not  true.  It  would  be  interesting 
if  one  could  estimate  the  diffusion  in  approximate  solutions  generated  by  finite 
difference  schemes  for  equations  in  classes  I  or  II. 

We  now  summarize  and  simplify  the  results  in  [17,18]  regarding  systems  of 
conservation  laws  with  coinciding  shock  and  rarefaction  curves. 

Let 

Ut  +  F ( U)x  =  0  (1) 

denote  an  n  x  n  system  of  conservation  laws.  Let  (X,R)  denote  a  real 
eigenfield  for  dF,  meaning  that 

dF  •  R  *  XR  , 

X  t  R.  Let  M  denote  any  1-dimensional  C3-manifold  in  U-space.  We  say  that 
the  shock  and  rarefaction  curves  coincide  on  M  if  M  is  an  integral  curve  of 
an  eigenvector  R,  and  for  every  Ui,U2  f  M,  is  in  the  Hugoniot  locus  of  U2 


for  some  scalar  s.  The  following  theorem  characterizes  the  coincidence  of  shock 
and  rarefaction  curves: 

THEOREM  1 .  The  following  are  equivalent  regarding  a  1-dimensional  C3 
manifold  M: 

(i)  The  shock  and  rarefaction  curves  coincide  on  M. 

(ii)  M  is  either  a  straight  line  in  U-space  or  else  simultaneously  an 

integral  of  R  and  a  level  curve  of  X  for  (X,R)  a  real  eigenfield 
for  dF . 

(iii)  System  (1)  reduces  to  a  scalar  conservation  law  on  M. 

Here  condition  (iii)  is  explained  by  the  following  definition: 

DEFINITION .  System  ( 1 )  reduces  to  a  scalar  conservation  law  on  a  1-manifold  M 
in  U-space  if  in  a  neighborhood  of  each  point  on  M  there  exists  a  coordinate 
system  <p  and  a  scalar  conservation  law 

ut  +  f(u)x  -  0  (2) 

such  that,  if  U  =  <p(u),  u(x,t)  is  a  weak  solution  of  (2)  if  and  only  if 
<p  •  u(x,t)  is  a  weak  solution  of  (1). 

For  condition  (iii),  f'(u)  *  X  •  ?(u)  in  the  case  that  the  coordinate 
system  u  is  taken  to  be  a  component  of  U  [18]. 

The  main  step  in  the  proof  of  Theorem  1  is  the  following  simple  lemma.  Let 
U(£),  <  £  <  ?2*  be  a  regular  CJ-parameterization  of  a  C  1-manifold  M 

such  that  U(£)  =  a(£)R,  a  f  R,  where  (X,R)  is  a  real  eigenfield  for  dF.  Let 
R(£)  =  R(U<£)),  X ( C )  =  X(U(£) ) . 

LEMMA.  If  U(£)  lies  in  the  Hugoniot  locus  of  a  point  UL  for  all 
£  e  then 

*(C) [U(5)  -  UL1  •  N(£)  -  0  (3) 

for  all  £  f  (£-*£,)»  and  all  vectors  N(£)  such  that  N(£)  •  R(£ )  ■  0. 


Thus  condition  (3)  says  that  either  X  ■  0  and  X  is  constant  on  M,  or 
[U(£)  -  UL]  •  N(£)  “  0  for  all  £  and  N(£),  which  implies  that  M  is  a 
straight  line  in  U-space  [18]. 

Proof .  Let  [U]  =  U(£ )  -  U^,  [F]  *  F(U{£))  -  F(UL).  Then  by  assumption 

s[U]  «  [F]  ,  (4) 

where  s  =  s(£)  is  a  function  of  £  since  [U]  and  [F]  are. 

Differentiating  (4)  gives 

s[U]  +  sU  -=  X&  ,  (5) 

and  differentiating  again  gives 

g[U]  +  2sU  +  sU  =  Xfi  +  XU  .  (6) 

Taking  the  inner  product  of  both  sides  of  (5)  with  N(£)  gives 

S[U(£)  -  n  ]  •  N(£)  -  0  .  (7) 

Xi 

Now  by  (7),  if  [u]  •  N(£)  ?  0  (i.e.,  [U]  is  not  parallel  to  &),  then 
s  =  0,  and  (5)  gives  s  =  X.  Putting  these  two  conditions  into  (6)  yields 

Xu  -  g[U]  . 

Thus  whenever  [U]  •  N(£)  ^  0,  we  must  have  X  *  0,  and  so  we  can  conclude  that 

X[U<£)  -  U_]  •  N(£)  -  0 

Li 

for  all  £. 

We  now  restrict  our  attention  to  2x2  systems  (1).  In  this  case,  let 
U  =  (u,v),  F  =  (f,g),  and  let  p,  q  denote  Riemann  invariants  for  (1),  by  which 
we  mean  any  C  function  on  U-space  whose  level  curves  are  integral  curves  of  a 
real  eigenvector  field  of  dF. 

DEFINITION .  System  ( 1 )  has  a  contact  field  in  a  neighborhood  U  of  U-space  if 
there  exists  a  Riemann  invariant  q  defined  in  U ,  Vq  f  0,  such  that  the  level 
curve  q  *  X  is  an  integral  curve  of  R  for  some  eigenfield  (X,R)  of  dF. 

Note  that  Vq  +  0  implies  that  any  system  having  a  contact  field  is  not 
linear,  and  is  not  a  nonlinear  transformation  of  a  linear  system. 


DEFINITION.  System  ( 1 )  has  a  line  field  if  there  exists  a  Riemann  invariant  q, 

(Jy 

Vq  f  0,  such  that  q  -  const,  is  a  line  of  slope  q  -  ■gjjf  i.e., 

qu  +  qqv  “  0  •  (8) 

Note  that  Vq  f  0  implies  that  any  system  having  a  line  field  does  not  uncouple 
in  this  field.  (In  an  uncoupled  field,  the  integral  curves  are  straight  lines  of 
constant  slope. ) 

We  prove  the  following  simple  theorems  which  give  necessary  and  sufficient 
conditions  for  a  2x2  system  to  have  a  line  or  a  contact  field. 
theorem  2.  A  2x2  system  of  conservation  laws  has  a  contact  field  if  and  only 
if 


f  -  uq  +  F(q) 

(9) 

g  ■  vq  +  G(q) 

(10) 

for  some 

smooth  functions  F 

and  G. 

PROOF. 

The  function  q  is 

a  Riemann  invariant  of  a  contact  field  if 

and  only  if 

qtu]  -  [f]  , 

(ID 

q[v]  -  [g]  , 

(12) 

for  all 

jumps  at  q  *  const. 

But  (11)  and  (12)  hold  if  and  only  if 

[f  -  qu]  *  0  , 

(13) 

[g  -  qv]  -  0  , 

(14) 

for  all  jumps  at  q  *  const.;  and  this  is  equivalent  to 

f  -  uq  +  F(q)  , 
g  *  vq  +  G(q)  . 

THEOREM  3.  A  2x2  system  of  conservation  laws  has  a  line  field  if  and  only  if 

g  *  fq  +  H(q) 


for  some  smooth  function  H,  and  some  smooth  solution  q  of  Burgers  equation 


PROOF.  A  smooth  solution  g  of  (15)  is  a  Riemann  invariant  if  and  only  if  the 


value  of  q  at  (u,v)  is  the  slope  of  the  integral  curve  of  a  line  field 
through  (u,v).  Thus  by  Theorem  1,  a  smooth  solution  g  of  (15)  is  a  Riemann 
invariant  for  system  (1)  if  and  only  if 

s[u]  =  [f]  , 
s[v]  =  [g)  , 


for  some  scalar  s,  whenever 


[v] 

tu]  ~  q 


But  this  condition  is  eguivalent  to 


q  -IaL 

q  [f] 


for  all  jumps  at  q  =  const./  which  holds  if  and  only  if 

[gf  -  g]  =  0 

for  all  jumps  at  q  =  const.  This  final  condition  is  eguivalent  to 

g  =  qf  +  H(q) 

for  some  smooth  function  H  and  some  smooth  solution  q  of  (15). 

The  following  corollaries  characterize  the  systems  in  classes  I,  II  and  III: 
COROLLARY  1 .  A  2x2  system  of  conservation  laws  has  two  independent  contact 
fields  with  Riemann  invariants  p  and  q  if  and  only  if 

f  =  uq  +  F^fq)  =  up  +  F2(p)  *  (16) 

g  =  vq  +  G^(q)  «  vp  +  G2(p)  ,  (17) 

for  some  smooth  functions  F^,  G^,  7p  ^  Vq.  Conditions  (16)  and  (17)  hold  if  and 
only  if  p  and  q  as  functions  of  u  and  v  are  given  by  the  inverse  of  the 
map 


/F^q)  -  F2(p)\ 

\G  (q)  -  G  (p )J 


(  18) 


PROOF.  Statement  (16)  and  (17)  follow  directly  from  Theorem  2,  and  (18)  is 


obtained  by  solving  for  u  and  v  in  (16)  and  (17). 

We  can  use  (18)  to  construct  a  2x2  nonlinear  coupled  system  of 


conservation  laws  having  two  contact  fields  as  follows:  choose  F^(w)  -  w2 


G^tw)  *  w3  for  i  *  1,2.  In  this  case 


Lv  W  ♦  pq  ♦  pV 


(19) 


We  can  invert  (19)  by  obtaining  p  =  -(q  +  u)  from  the  first  equation  in  (19) 
and  substituting  this  into  the  second  equation.  By  solving  the  resulting 
quadratic  equation  in  q,  we  obtain 


(P) 

V 


»(: 


4v  -  3u‘ 


A 


4v  -  3u 


(20) 


Choosing  -  for  p  and  +  for  q  in  (20),  we  obtain  f  =  -u  -  v  and 


g  =  u3  +  uv  by  substituting  (20)  into  (16),  (17).  Thus  the  2x2  system 


u..  -  (u2  +  v). 


(21) 


”  0  , 

Vt  +  (u3  +  uv)x  =  0  , 
has  two  contact  fields  with  wave  speeds  p  and  q  given  by  (20). 

COROLLARY  2 .  A  2x2  system  of  conservation  laws  has  one  contact  field  and  one 
line  field  if  and  only  if 

f  =  (u  +  a)q  ,  (22) 

g  =  (v  +  b)q  ,  (23) 

for  some  constants  a,  b,  and  some  smooth  function  q  of  u  and  v.  In  this 


case  q  is  the  wave  speed  and  Riemann  invariant  for  the  contact  field,  and 

v  +  b 


P  = 


u  +  a 

is  the  Riemann  invariant  of  the  line  field  which  satisfies  (15). 


-7- 


(24) 


PROOF .  (See  (2. 1 0 )- ( 2 . 18)  of  [17].) 


Systems  of  type  (22),  (23)  with  a  =  b  =  0  have  arisen  in  applications,  and 
were  studied  in  [5,8,16]. 

COROLLARY  3.  A  2x2  system  of  conservation  laws  has  two  independent  line 
fields  if  and  only  if 


f  = 


H,(p)  -  H2(q) 


q  -  p 


g  = 


qH^p)  -  pH2(q) 


q  -  p 


(25) 


(26) 


for  some  smooth  functions  H-j  and  H2 ,  and  two  smooth  independent  solutions 
p  and  q  of  (15). 

PROOF .  By  Theorem  3,  system  (1)  has  two  independent  line  fields  if  and  only  if 

g  =  fp  +  H-j  (p)  (27) 

and 


g  =  fq  +  H2(q)  (28) 
for  smooth  independent  solutions  p  and  q  of  (15).  Thus  solving  for  f  and 
g  in  (27),  (28)  gives  (25),  (26). 


The  system 


«t  ♦  {  T+TTT  1 

X 

vt  -  <  } 

X 


0  , 


0  , 


for  0  <  k  <  1,  u  >  0,  v  >  0,  is  an  example  of  a  system  in  class  III  arising  in 
the  study  of  multicomponent  chromatography  [1,14]. 

Acknowledgement:  Thanks  to  Rich  Sanders  for  working  out  example  (21)  with  me. 
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